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It is shown that the T, operator used in a previous study to define the open-shell 
D, diagnostic is invalid, and leads to an arbitrary definition of the open-shell D, 
diagnostic. A new definition is proposed that eliminates this ambiguity and 
approximately restores the mathematical relationship previously noted between the 
closed-shell D t and T, diagnostics. Statistical comparison of the T, and D, diagnostics 
shows a very high degree of correlation between them, although it is argued that both 
diagnostics used together can provide more information than either can separately. 


Recently, Leininger et al. [1] defined the open-shell version of the D, diagnostic, 
the closed-shell version of which was defined by Janssen and Nielsen in 1998 [2] The D, 
diagnostic is based on the matrix 2-norm of the coupled-cluster t, amplitudes, and, like the 
T, diagnostic [3-6] which is based on the Frobenius norm, it is designed to give an 
indication of the quality of results to be expected from a singles and doubles coupled- 
cluster (CCSD) calculation. Janssen and Nielsen showed that there was a rigorous 
mathematical relationship between the D, and T, diagnostics for closed-shell theory, 
which is due to the fact that they are both based on the coupled-cluster ti amplitudes. 

For the open-shell version, however, no such relationship exists even though both die D, 
and T, open-shell diagnostics were defined for coupled-cluster wavefunctions based on 
restricted Hartree-Fock reference functions. There is a difference, however, since the 
open-shell Ti diagnostic was defined using an open-shell coupled-cluster wavefunction 
that is an S x eigenfunction [6] based on symmetric spin-orbitals [7], whereas Leininger et 
al. defined a new T/ operator based on symmetric and antisymmetric combinations of 
annihilation and creation operators. This T, operator was not actually used in calculating 
the wavefunction, rather the standard spin-orbital equations were used, and the 
resulting t, amplitudes transformed into those for the newly defined T, operator. We 
show below, however, that the new T, operator is not valid and that the relationship 
between the S z open-shell t, amplitudes and the new amplitudes is arbitrary. Hence, the 

definition of the open-shell Di diagnostic is also arbitrary. 

One of the purposes of this study is to eliminate the ambiguity in the open-shell 
Dj diagnostic by basing it on the restricted open-shell coupled-cluster wavefunction 
developed in Ref. [6]. The mathematical relationship between the D, and T, diagnostics 
is now valid for both the open- and closed-shell versions. We provide test examples of 


1 tjlee@mail.arc.nasa.gov 


Relationships between the coupled-cluster amplitudes in the S, and the usual & basts 
were also given in Ref. [6], For current purposes, a brief summary ts given here. / 

operator, spin rotated into an S : basis, is given by Eq. (45) m Ref. [6] and is 

r, = f-E m + v"A ai + 4l /; [(aa) + (xa)] + 42 f* [(x0)'(W)] , (3) 

where f and v are the t, amplitudes in the S, basis. The f amplitudes correspond to orbital 
relaxation parameters whereas the v amplitudes correspond to a type o _doub e ■excitoti 
in the basis (see Ref. [6] for a more complete discussion). Equating Eq. (1) and q. ( ) 
allows us to determine relationships between the f and t amplitudes: 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


fr=(c+*S)i 2 • 

r x =ci4i , 

f*=t*$l42 , 

v ?=( C -'?)' 2 . 

where Eq. (7) is given for completeness. Hie open-shell D, diagnostic is *en defined* 
the matrix 2-norm of the f amplitudes analogous to the procedure used in Ref. [1]. That 
is the three open-shell amplitude blocks are treated separately and the maximum matnx 
2-nomi is taken as die D, diagnostic in order to preserve die energy invar, ance properties 
of the open-shell D, diagnostic (see Ref. [1]). 

D 1 =max(||y; B || 2 ,||/;|| 2 ,|y:i 2 ) • (8) 

Again, following Ref. [1], the matrix 2-norm is computed according to 

ML = ” ^max 5 

where X is the largest eigenvalue of the matrix RR T and <r„ is the largest singular 
value of R. Similar to the relationship between the open- and closed-shell ^‘agnostics, 
if there are no open-shell orbitals *, then the open-shell D, diagnostic is identical to the 
closed-shell definition [2], Defining the open-shell D, diagnostic in tins way has the 
additional advantage that the equation relating the closed-shell D, and T, diagnostics, 

■J2T < D now approximately holds for the open-shell diagnostics in e inn 
number of open-shell electrons being much smaller than the total number of elects 
Test case examples of the new D, diagnostic are given in the next section followed by 
critical comparison to the 7) diagnostic. 
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Results and Discussion 


A. 


New Open-Shell Di Values 

Test case examples of the new open- 


shell D, diagnostic together with the open- 



both the closed-shell and open-shell molecules included in these two studies To give an 
idea of how strong the correlation between T ■ and D; is, we can compare o e 
correlation coefficient between the T, and * diagnostics is the P— e 
equivalent of the T, diagnostic [1 1]), as well as to the correction * e 

coupled-cluster and perturbation theory D, diagnostics [2]. Using between 

shell molecules contained in Ref. [2], we obtain a correlation coefficient of 0.98 betwe 
the and 5, values, and a correlation coefficient of 0.95 between the < 
perturbation foeoty D, diagnostics. This shows that the correlati, an between i te A an 
T, coupled-cluster diagnostics is almost as strong as the correlation with their respect 

^"^"C^Son is to be expected because of dte similarities 
between the L diagnostics. Both the T, and D, diagnostics are based on the orbital 
relaxation parts of the t, amplitudes, both were designed to exhibit the same orbital 
invarianceproperties that the CCSD energy possesses, and in spite of su «8« s ““ “ 
contrary n,2],hey were both designed to exhibit the mathematical property of size- 
“Si For example, bod, diagnostics will yield the same value for a single helium 
"for any number rf non-interacting helium atoms. However, it was pointed out m 
Ref. [5] that the T, diagnostic may fail to indicate that a small region o a 
difficult to describe properly if the rest of the molecule is well described ' * * e C “ D 
level of theory and in Ref. [2] Janssen and Nielsen gave a numerical examp e of this 
pta— In essence, L t, diagnostic is an average over the whole molecule and foe 

contribution from the small problem area is swamped by that | 0m ' f^^mostic or 
molecule which is well described. This can be viewed as a failure of the T, diagnostic or 
it can ato be viewed as a success since the majority of the molecule is well described and 
the T, diagnostic indicates this. Conversely, the D, diagnostic is designe to yie a arge 
value for ftage molecule with only a small problem area. Again, this could be viewed as 
™s orTkilure since most of die molecule is well described, which the Diagnostic 
does not suggest, but there is one problem area which the D, diagnostic does indica e. 
Another simation where the two diagnostics could give conflicting information is when 
^"cumulation of correlation effect wherein there is no one orbrnt retotion 
parameter or excited state that is very important, but rather there are several states or 
orbital relaxation parameters that are moderately important. In tins case the T, diag 
would be larger than usual, indicating the need to treat higher-order correlation wlu 
would be missed by the D, diagnostic. The important point, is t^t togethertte^ m ^ 
Di diagnostics provide more information than either does one, 

We stress, however, that we also believe that diagnostics based on two ‘l >a " lde 

components ofthewavefunction, such as the D, diagnostic [12], can provi ea l ion , 

important information in assessing the quality of a particular calculauon 

Another quantity to consider is the ratio T,ID„ and these are included m Table 1 
for the open-shell molecules studied here. In a perfectly homogeneous ^^m.tiusratio 
can be derived from the approximate mathematical relationship between the T, and D, 
diagnostics to be l/S ■ A homogeneous system in this se^e is def ‘“ d “ 
foe contribution from all of the f amplitudes is identical. Thus, the further the T,ID, ratio 
deviates from l/V2 , the greater the non-homogeneity of the electronic structure of the 
molecule being studied. We note that the T,iD, ratio will deviate from a perfectly 
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